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The calculation of higher twist (or dimension) corrections to physical quantities using operator product ex- 
pansions is delicate. If dimensional reguiarization is used to regulate the ultra-violet divergences then there are 
ambiguities in the Wilson coefficient functions due to infra-red renormalon singularities. With a hard ultra violet 
cut-off, such as the inverse lattice spacing a, there are no renormalon ambiguities, as a result of cancellations be- 
tween terms which in finite orders of perturbation theory diverge as inverse powers of a, and those which diverge at 
most logarithmically. In this lecture I review these questions, explaining the steps necessary to obtain predictions 
for physical quantities from lattice measurements of matrix elements of higher dimensional operators. The ideas 
are illustrated by considering quantities computed using the heavy quark effective theory beyond leading order 
in the heavy quark mass. 



1. Introduction 

In recent months there has been renewed inter- 
est in the evaluation of power corrections to QCD 
predictions for hard scattering processes and re- 
lated quantites, and in particular in attempting 
to understand the implications of the presence 
of renormalon singularities. Among the impor- 
tant physical quantities which are being studied 
are the higher-twist contributions to deep inelas- 
tic structure functions, the l/mg corrections to 
masses and operator matrix elements in heavy 
quark physics (where mg is the mass of the heavy 
quark Q), and corrections to predictions for event 
shape variables in jet physics and to the cross- 
section for the Drell-Yan process (i.e. the produc- 
tion of lepton pairs with a large invariant mass in 
hadronic collisions). For some of these processes 
the discussion can be formulated in terms of the 
operator product expansion, and it is such pro- 
cesses which will be considered in detail in this 
talk. The relevance of renormalons in field the- 
ory, and for operator product expansions in par- 
ticular, has been developed in the papers listed 
in ref.[l]. This work has been exploited and ex- 
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tended to the study of renormalons in the Heavy 
Quark Effective Theory (HQET) in refs.[2]-[6]. 
There has however also been some work recently 
on the role of renormalons in hard scattering pro- 
cesses for which there is no operator product ex- 
pansion (see refs [7]-[12] for example). In this talk 
I will discuss the implications of the existence of 
renormalons in general, and in evaluating physi- 
cal quantities from lattice simulations in particu- 
lar. This may seem a surprising subject for a lec- 
ture, since renormalons are a consequence of the 
divergent behaviour of perturbation theory, and 
lattice simulations are a non-perturbative tech- 
nique. It should be remembered however that in 
determining physical quantites from operator ma- 
trix elements computed on a lattice, one has to go 
through the process of "matching" which is per- 
formed in perturbation theory. It is at this stage 
that the question of the existence and importance 
of renormalons arises. The study of renormalons 
in lattice field theory, with its "hard" ultaviolet 
cut-off a^^ (where a is the lattice spacing), also 
provides important general insights, complement- 
ing the usual studies which are almost exclusively 
carried out using dimensional reguiarization of 
ultra-violet divergences. 

In lattice field theory higher dimensional oper- 
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ators mix with lower dimensional ones with the 
same quantum numbers, with mixing coefficients 
which diverge as inverse powers of the lattice 
spacing. Thus for example, the kinetic energy 
operator in the HQET, hD^h (where h represents 
the field off the heavy quark) mixes with l/a^hh 
and l/ahDih. These power divergences must be 
subtracted non-perturbatively in order to define 
"physical" matrix elements, since factors of the 
form 



exp 
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QCD 
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combined with inverse powers of the lattice spac- 
ing can give finite contributions [13]. Renor- 
malons are an example of such non-perturbative 
effects. Renormalons are also present in calcula- 
tions using dimensional regularization for those 
quantities for which power divergences would be 
expected by power counting. Using dimensional 
regularization however, power divergences and 
the perturbative mixing of operators of different 
dimensions are absent because of the lack of a 
hard cut-off. 

For most of this lecture I will consider compu- 
tations in the HQET, however the principal fea- 
tures concerning the appearance and cancellation 
of renormalons are more general. Thus the dis- 
cussion below can readily be generalized to other 
applications of operator product expansions. In 
the HQET one computes physical quantities as 
series in inverse powers of the heavy quark mass. 
The Dirac term in the QCD action is replaced by 

Q{i^ - mQ)Q hy{iv ■ D)K + ^-^KiiOfK 

+ /i«<T„^J^"^/i. + 0{l/ml) , (2) 

where Q and hy represent the fields of the heavy 
quark in QCD and the HQET respectively, and 
V is the quark's four-velocity. In all the simu- 
lations described below we will take the quark 
to be at rest, and will denote the correspond- 
ing heavy quark field by h (without any sub- 
script). Cmag is a constant determined by match- 
ing the effective theory onto the full one (QCD). 
The corresponding constant for the kinetic term 



{hy{iD)'^hv) is equal to one by reparametrization 
invariance [14]. In addition to the expansion of 
the action in eq.(2), one must also expand the 
operators whose matrix elements are being com- 
puted (this will be discussed below). For a com- 
prehensive review of the structure and applica- 
tions of the HQET, and references to the original 
literature, see ref.[15], and for discussions of lead- 
ing order computations using lattice simulations 
of the HQET see refs.[16,17]. 

Much of the discussion in this lecture will be 
based on the two interesting papers which pointed 
out the presence of renormalons in the HQET, 
and stressed that the pole-mass of a quark cannot 
be a physical parameter and should not be used as 
the expansion parameter [2,3]. The applications 
to the lattice formulation of the HQET, and the 
numerical results presented below, are based on 
work carried out with Guido Martinelli, Marco 
Crisafulli and Vicente Gimenez [18,19]. 

The plan of this lecture is as follows. In the 
next section I present the main results concern- 
ing the appearance of renormalons in operator 
product expansions and their significance, delay- 
ing the corresponding detailed explanations until 
the later sections. In sections 3 and 4 I discuss the 
appearance of renormalons in Borel transforms 
of perturbation series, and present calculations 
of their effects in the limit of a large number of 
light quark flavours. The presence of renormalons 
in the pole mass in demonstrated in section 5. 
Section 6 contains a discussion of renormalons in 
the HQET regulated with a hard ultra-violet cut- 
off, and in the lattice formulation of the HQET 
in particular. Although the pole mass is not a 
physical parameter, short-distance definitions of 
the quark mass are physical, and in section 7 I 
explain how these can be evaluated from simu- 
lations in the HQET, and present some results 
for the MS mass. The evaluation of the kinetic 
energy of the heavy quark is explained in sec- 
tion 8. This calculation demonstrates the general 
non-perturbative techniques needed to subtract 
power divergences and renormalon uncertainties 
from matrix elements of higher dimensional op- 
erators which were proposed in refs.[18,19]. In 
section 9 I discuss a possible definition of the pa- 
rameter A, the difference of the hadron and heavy 
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quark masses, which is commonly introduced in 
applications of the HQET. Numerical results for 
A are also presented. Finally section 10 contains 
a summary and my conclusions. 

2. Overview 

In this section I will present the principal re- 
sults concerning the appearance and cancellation 
of renormalon singularities in operator product 
expansions, using the Heavy Quark Expansion as 
an illustration. Consider the evaluation of the 
matrix element of some local QCD operator con- 
taining one or more heavy quark fields. We call 
this operator 0*5CD p^j, example if we wish 
to evaluate the leptonic decay constant /a, then 
qQCB Tj^Q^if^ i-,g ^jjg axial current b'y^^'y^q, where q 
represents the field of the light quark. Using the 
HQET we expand 0*°^*^^ in inverse powers of the 
mass of the heavy quark itiq 
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C2(mQ///,)0^Q=T(M) + 0(l/m^) (3) 



where fi is the scale at which the operators of 
the HQET, O^^^^ and O^^^'^ , are renormal- 
ized. We consider here the simple situation for 
which there is a single operator in each of the 
first two orders of the expansion, but the discus- 
sion below can be easily generalized to the case 
in which there are more operators. 

When evaluating matrix elements of qQ^id 
yond the leading order in the 1 /mg expansion, in 
addition to the corresponding higher dimension 
operators on the r.h.s. of eq.(3), it is also nec- 
essary to keep higher order terms in the heavy 
quark action, as in eq.(2). Throughout the dis- 
cussion below it is implied that this is done. 

In general the QCD operator O'^'^'^ also re- 
quires renormalization and is defined at some 
scale M. Unless specificahy needed I will sup- 
press the dependence on M in O^'^^ and in the 
coefficient functions C^. 

I will now summarize the main points which 
I wish to make in this lecture, postponing a de- 
tailed discussion to the following sections: 

• If we restrict the calculation to the lead- 
ing order in l/mg, i.e. if we neglect the 



0(1 /niQ) terms, then the perturbation se- 
ries for the coefficient function Ci diverges, 
and is not Borel summable. As wih be ex- 
plained below this is due to a singularity in 
the Borel transform of the perturbation se- 
ries, called an infra-red renormalon. Thus 
there is an ambiguity in the evaluation of 
Ci , coming from the different possible ways 
of defining the series. This ambiguity is of 
Oil/ruQ). 

• Formally therefore, we should not include 
the 0{l/mQ) corrections in eq.(3) until we 
have computed sufficiently many terms in 
the perturbation series for Ci to control its 
divergent behaviour. 

• For this talk I will restrict the discussion 
to problems for which the matrix elements 
of oQ'^° have no renormalon ambiguities. 
In general such non-perturbative effects do 
exist, and appear on the right hand side of 
eq.(3) in the matrix elements of the opera- 
tors of the HQET. They do not affect the 
coefficient functions. 

• In renormalization schemes based on the 
dimensional regularization of ultra-violet 
divergences in the HQET, such as the 
MS renormalization scheme, the matrix 
elements of O^^^^ are also not Borel 
summable in perturbation theory, due to an 
ultra-violet renormalon singularity in their 
Borel transforms. The corresponding am- 
biguities in the matrix elements of O^^^^ 
cancel those in Ci. 

• If a hard iiltra-violct cut-off is used, such 
as the lattice spacing in the lattice HQET, 
then the matrix elements of O^*^^^ do not 
have ambiguities due to ultra-violet renor- 
malons. This, in turn, imphes that the co- 
efficient functions Ci do not contain ambi- 
guities due to infra-red rcnormalons. In the 
lattice theory it is natural to present the dis- 
cussion in terms of bare operators in the ef- 
fective theory, defined with the lattice spac- 
ing a as the ultra-violet cut-off. Through- 
out this talk I will assume that mqa 3> 1. 



Of course if the inverse lattice spacing was 
much smaller than the heavy quark mass 
then there would be no need to use the 
HQET. 

The absence of renormalon ambiguities in 
Ci in the lattice theory is due to a cancella- 
tion between leading and non-leading terms 
in perturbation theory, i.e. terms which be- 
have logarithmically with mqa and non- 
leading ones of 0{l/mQa). 

In the lattice theory, matrix elements 
of higher dimensional operators (such as 
qHQEt-j jjypj,gg gg inverse powers of the lat- 
tice spacing, and are hence manifestly un- 
physical. For example 

{H\hD^h\H) 0{lla^) (4) 

where H represents the heavy hadron. 

The subtraction of the power divergences in 
perturbation theory introduces renormalon 
ambiguities. Thus for example, the per- 
turbation series of the terms which diverge 
quadratically, i.e. those which are 0{l/a^), 
in eq.(4) contains a renormalon ambiguity. 

The perturbation series for the pole mass 
has a renormalon ambiguity of O(Aqcd) [2, 
3]. This is not the case for short-distance 
definitions of the heavy quark mass, such 
as 

ruQ^mfimf). (5) 

It is possible to compute mg (and other 
short-distance definitions of the mass), us- 
ing only simulations in the HQET. In 
ref.[19] it was found that 

rUb = 4.17 ± 0.05 ± 0.03 GeV -I- 0(l/mb) .(6) 

The presence of rcnormalons in physical 
quantities for which there is no Operator 
Product Expansion, such as the Drell-Yan 
process or in event shape variables in jet 
physics, is a subject currently under inten- 
sive investigation [7]-[12]. It is hoped that 



these studies will provide important phe- 
nomenological information about the sub- 
asymptotic (non-leading twist) behaviour of 
physical quantities. In lattice QCD there 
are analogous questions about the pres- 
ence of ambiguities in perturbation series 
for quantities which contain power diver- 
gences, e.g. is there an ambiguity in the 
perturbative evaluation of the critical mass 
when using Wilson fermions? I will return 
to these questions at the end of this talk. 

3. The Borel Transform and Renormalons 

In this section I define the Borel transform of a 
perturbation series and explain what is meant by 
a renormalon singularity. Consider the perturba- 
tion series for some quantity F: 

Fia,) = Y,Fn[^aM) (7) 

where, for later convenience, we choose to in- 
clude /3o/47r in the expansion parameter (/Jq = 
11 — 2/3n/ is the first term in the perturbative 
expansion of the /3-function) . The Borel transorm 
F{u) of F{as) is defined by 

F(u) = Fo6(u) + V -F^+iU^ . (8) 

n=0 



If the theory is Borel summable then F is recon- 
structed by taking the Laplace transform 




F(as) is the unique function which satisfies cer- 
tain analyticity conditions in the Qj-planc with 
the asymptotic expansion (7) (see ref.[20] for ex- 
ample) . 

If the coefficients Fn grow like nl or faster, then 
the factor of n! in the denominator of the terms on 
the r.h.s. of eq.(8) will not be sufficient to prevent 
singularities appearing in F{u) for values of u on 
the positive real axis. In such cases the integral 
in eq.(9) is not defined. Of course the contour of 
integration in eq.(9) can be deformed away from 
the singularities, but in general the resulting inte- 
gral will depend on how this is done, leading to an 
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Figure 1. One loop graph contributing to the 
total cross section for e+e~ annihilation. The 
dashed line repesents a virtual photon or Z- 
Boson. 



ambiguity in the result. Renormalons is the name 
given to the singularities on the positive real axis 
in the u-plane. 

In order to obtain some intuition concerning 
the appearance of renormalons consider the one 
loop graph in Fig.l, which can be considered, for 
example, contribution to the amplitude for 
e+e^ hadrons. It may be expected that part 
of the effect of higher loop corrections will be to 
renormalize the coupling constant, and so we take 
the coupHng constants at the vertices to be g{k), 
where k is the momentum in the gluon propaga- 
tor. Now expanding g{k) as a power series in 
where /j is some fixed renormalization scale, we 
readily find that the coefficients in perturbation 
series grow Hke n\. 

In general the renormalon singularities are 
branch points of cut singularities, whose position 
and nature are determined by the renormaliza- 
tion group, but not their residues. In some sim- 
plified examples below, motivated by the large 
Nf expansion, where Nf is the number of quark 
fiavours, they appear as simple poles: 



F{u) = y + 



U — Ui 



(10) 



where we choose the ordering Q < Ui < U2 < 
U3 ■ ■ ■ . Writing 



u-Ui (u-ui).p 



+ r]iS{u-Ui) 



(11) 



where the subscript V stands for "principal part" , 
and rii defines the prescription used to define the 
pole. We define AF, the ambiguity in F as the 
coefficient of rji [5] , 



AF = n exp 



—Anui 
/3oas{id) 



ri 



QCD 



.(12) 



From eq.(12) we see that the ambiguity is indeed 
"exponentially small" , and hence is suppressed by 
a power of Aqcd/a*- 

4. The Large Nf Limit 

Of course in general it is not possible to sum 
all higher order graphs completely. Much insight 
can be obtained by studying the behaviour of per- 
turbation series in the limit of a largo number of 
fight quark fiavours Nf. In this fimit the gluon 
propagator is given by the sum of bubble dia- 
grams shown in Fig. 2, since it is only for these 
graphs that each additional factor of ctg is accom- 
panied by a factor of Nf. The geometric series for 
the gluon propagator in this limit can readily be 
summed, leading to the following expression for 
its Borel Transform in the Landau gauge: 



k^'k" - k'^g^"' 



(13) 



where a and 6 are colour labels. C is the finite 
(renormalization-scheme dependent) term in the 
one-loop graph, which is given by 



-^{H-kV,') + c) 



(14) 



In the large Nf limit asymptotic freedom is 
lost, and instead we replace the factor — 2/37V/ 
in eq.(14) by /3o. In this way we sum the leading 
powers of iSq in each order of perturbation theory, 
and the procedure is motivated by the intuition 
that it is the scale dependence of the coupling 
constant which gives rise to renormalons, at least 
partially. 

In the approximation described above, all the 
dependence on the coupling constant comes from 
diagrams with a single resummed gluon propaga- 
tor. Thus the 0{1/Nf) contribution to the Borel 
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Figure 2. The sum of diagrams contributing to the gluon propagator in the large Nf limit. The "bubbles" 
represent loops of light quarks. 



tion for the pole mass, however, as 




p + k 



Figure 3. Diagram contributing to the leading 
non-trivial term in the Borel transform of the 
quark propagator in QCD. The box containing 
u implies that for the gluon propagator we take 
the expression in eq.(13). 



transform of the quark propagator is given by the 
diagram in Fig. 3 (this is the leading non-trivial 
contribution). In the following section wc study 
the graph in Fig. 3, and demonstrate the appear- 
ance of renormalons in the Borel transform of the 
"pole-mass" . 



5. Renormalons in the Pole Mass 

Wc are now readjf to study the Borel transform 
of the perturbation series for the pole mass. The 
pole mass is given by 



m 



pole = niQ + I]( 



(15) 



where mg is the bare heavy-quark mass and S is 
the self-energy diagram, whose Borel transform is 
represented in Fig. 3. Eq.(15) is an implicit equa- 



mpoie =mQ + yj^j ' (16) 

in the large Nf limit we can replace mpoic with 
TTiQ on the right hand side of eq.(15). Thus the 
Borel transform of the pole mass can be obtained 
by simply evaluating the graph in Fig. 3, on shell, 
i.e. with = tUq. 

It is straightforward to see that on-shell there 
is a pole singularity at u = 1/2. The infra-red 
behaviour of the graph in Fig. 3 is given by 

/ ^'^715^ (17) 

Jk- 



(fc2)l+-fc 

which, using dimensional regularization, has pole 
singularities of the form r(l — 2u), i.e. the first i.r. 
renormalon appears at u = 1/2. The correspond- 
ing ambiguity in the pole mass is of O(Aqcd) (see 
eq.(12) above): 

Ampoie = -^e-^/2^c 



l3o 



VQCD 



(18) 



For the following discussion it will be useful to 
observe that off-shell, i.e. for ^ mg there is 
no pole at u = 1/2. In this case the infra-red 
behaviour of the graph in Fig. 3 is given by 

1 



d'^k- 



(19) 



which is manifestly finite for u = 1/2. 

We are now ready to demonstrate the cancella- 
tion of renormalons described in section 2. Con- 
sider the expansion of the inverse heavy quark 
propagator in QCD, as a series in inverse powers 
of the heavy-quark mass. In perturbation theory 
it takes the form: 

Sp^ (p) = m- mpoie + C{mQ/ij)S^g {v ■ k, fi) 

+ 0{l/mQ) (20) 
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— ^ 




k + q 



Figure 4. Diagram contributing to the leading 
non-trivial order term in the Borel transform of 
the quark propagator in the HQET (the bold line 
represents the quark propagator in the effective 
theory). The box containing u implies that for 
the gluon propagator we take the expression in 
eq.(13). 



where S^g is the inverse of the heavy quark prop- 
agator in the HQET. The subscript P impHes the 
projection: 



(21) 



2 2 

k is defined hy p = mv + A;, where m will be 
chosen conveniently later. The freedom to choose 
m in this decomposition of p is equivalent to the 
freedom to choose a renormalization condition in 
the HQET to fix the mass counterterm. 

The Borel transform of eq.(20) takes the form: 

Sp^{p,u) = m{u) - mpoie(M) + 

C{mQlii, u) * S;s\v ■ k, fi, u) + Oil/niQ) (22) 

where the tilde represents the Borel transform 
and * the convolution. We have seen above that 
TTipoie has an i.r. renormalon at u = 1/2. S~g has 
an u.v. renormalon at u=l/2. This can easily be 
deduced by power counting. The graph in Fig. 4 
has the ultra-violet behaviour 

/ ^'?7^^ (23) 

./large q (Q ) Q 

which yields a pole at u =1/2. 

The i.r. renormalon in mpoie is cancelled by the 
u.v. renormalon in . This can be verified by 



explicit calculation, but follows directly from the 
observation made above, that the Borel transform 
of S~^{p) has no renormalon at m = 1/2 for off- 
shell values of the momenta p. 

The connection between the presence of power 
divergences in perturbation theory and renor- 
malons can also be demonstrated using power 
counting. The integral in eq.(23) has a logarith- 
mic ultra-violet divergence at u = 1/2, which im- 
plies that it has a power (in this case linear) di- 
vergence at u = 0. The one- loop contribution to 
the quark self energy corresponds to m = 0. Con- 
versely, if a low order graph with a single gluon 
propagator has a power divergence in perturba- 
tion theory, then the Borel Transform of the cor- 
responding scries of graphs with bubble insertions 
on the gluon propagator has singularities at pos- 
itive values of u, i.e. renormalon singularities. 

The appearance and cancellation of renor- 
malons in eq.(22) is an example of the general 
structure outlined in section 2. In this case we 
can view OQ^^ as the correspondinj|_ Dirac op- 
erator in QCD, and of^^^ and ^ as 
and hvV ■ Dh^ respectively. 

6. Renormalons and Lattice Perturbation 
Theory 

The discussion in the previous section was 
based on the use of dimensional regularization of 
the ultra-violet divergences. With a hard ultra- 
violet cut-off, such as the lattice spacing, there 
cannot be a renormalon singularity in S~g (see 
eq.(20) ). This is because the cut-off itself en- 
sures that the integral corresponding to eq.(23) is 
finite at all values of u. 

In lattice perturbation theory, the inverse prop- 
agator in the HQET takes the form 

S;^\v ■k)=vk[l + (_^^ \n{-2av ■ k) 



+d)] - as{a)— + 



(24) 



where jh is the one-loop contribution to the 
anomalous dimension of the heavy quark field, 
as{a) is the bare coupling constant and d and X 
are constants. The presence of the linearly di- 
vergent mass term, —as{a)^, is a manifestation 
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of the mixing of the operator h^v ■ Dhy with the 
lower dimensional one hyh^. 

As mentioned above the Borel transform of the 
scries in eq.(24) does not have any renormalon 
singularities. However the absence of these sin- 
gularities occurs through an interesting cancel- 
lation. The series of terms in eq.(24) which di- 
verge linearly does have a renormalon at m = 1/2. 
This renormalon cancels against that present in 
the Borel transform of the series of terms propor- 
tional to V ■ k. To see this consider the graph in 
Fig. 4 for V - k = Q. The infra-red behaviour of the 
graph is given by 



L 



(25) 



which has pole at u = 1/2, demonstrating the ex- 
istence of the renormalon in the series of terms 
which diverge linearly in each order of perturba- 
tion theory. Since at v ■ k ^ there is no singu- 
larity, the renormalon is cancelled by that in the 
terms which, in each order, are proportional to 
V ■ k. The structure of the Borel transform of the 
inverse propagator as u 1/2 is: 



Sj {v ■ k, u) 



1 - 2u 



V ■ k{av ■ k) 



-2u 



a 



(26) 



plus terms which are manifestly non-singular at 
u = 1/2. The residue of the pole at m = 1/2 
is zero as a result of the cancellation described 
above. 

7. The MS Mass 

Although the pole mass is ambiguous due to the 
presence of renormalon singularities, it is possible 
to evaluate any "short-distance" mass using lat- 
tice simulations in the HQET and perturbation 
theory. In this section we present a determina- 
tion of the MS mass mg defined in eq.(5). Con- 
sider the correlation function of some heavy-light 
hadron, for example that of the pseudoscalar me- 
son P, 

Cit) = X^(0|/i(x)7^g(x)g(0)7^ft(0)|0) (27) 
= exp(-ei) . (28) 



The matrix element Z is proportional to the de- 
cay constant fp in the static approximation, and 
has been the subject of many lattice computa- 
tions [17]. The quantity e is linearly divergent as 
a ^ 0. Now rriQ can be obtained from the com- 
puted value of e and the known value of the mass 
of the meson P (in practice that of the B-meson) 
by the perturbative relation [19] 



mp 



CtsX 



4c^ 



-I- 



(29) 



The perturbation series in each of the two fac- 
tors in eq.(29) have renormalon singularities at 
u = 1/2, and these singularities cancel in the 
product. The first factor corresponds, in pertur- 
bation theory, to the determination of mpoie from 
the lattice simulation of the effective theory. The 
second factor is the perturbation series relating 
mpoie to rfiQ . In this way mg can be determined 
from the computation of e. The remaining am- 
biguities in eq.(29) correspond to renormalons at 
u = 1, and hence are of 0(AQQj-,/mQ). Since in 
the present discussion we do not include l/mg 
corrections to the HQET action, these remain- 
ing ambiguities are of the same order as other 
uncertainties. In principle we can also systemati- 
cally reduce the ambiguities and uncertainties to 
higher orders in l/mg. 

We have recently determined Wtt, in this way, 
using the precise results for e obtained by the 
APE collaboration and found [19]: 



mb = 4.22(7) GeV at /3 = 6.0 
nib = 4.15(8) GeV at /3 = 6.2 . 



(30) 
(31) 



Combining the results in eqs.(30) and (31) we ob- 
tain our best estimate which has been given in 
eq.(6) above. 

In obtaining these results we have used one- 
loop perturbation theory in the matching equa- 
tion (29). It will be very important to extend this 
to higher orders to be confident that we control 
the cancellation of the renormalon singularities 
with sufficient precision. 
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8. The Kinetic Energy Operator 

In this section I discuss the evaluation of the 
matrix elements of the kinetic energy operator 
which is proportional to hD^h. These matrix ele- 
ments arc needed, for example, in plicnomeno- 
logical applications of the HQET to the spec- 
troscopy and inclusive decays of heavy hadrons. 
This calculation provided an illustration of the 
general techniques we propose for the subtrac- 
tion of power divergences, which appear due to 
the mixing of higher dimensional operators with 
lower dimensional ones with the same quantum 
numbers. 

Wc start with the observation that the matrix 
elements of the bare lattice operator are huge, e.g. 
from a simulation at /3 = 6.0 on a 16^ x 32 lattice 
using the SW-Clover quark action at k = 0.1425 
we found [19], 

a^ffi^ = -0.72 ±0.14, (32) 
2Mb 

which corresponds in physical units to 

ffi|».-2.9GeV^ (33) 
2Mb 

whereas for a physical matrix element we would 
expect the result to be of O(Aqqq). Of course 
the matrix element of the bare operator is not 
a physical quantity, indeed, as explained above, 
it is quadratically divergent as a 0. How- 
ever the result above illustrates that the neces- 
sary subtractions will be very large, and hence 
the determination of physical effects will be dif- 
ficult. In one-loop perturbation theory the term 
to be subtracted is — 5.19Qs/a^ which is of the 
right order of magnitude. However, the series of 
quadratically (and linearly) divergent terms have 
renormalon ambiguities and hence a subtraction 
of these divergences in perturbation theory does 
not lead to a physical matrix element. It is there- 
fore necessary to perform these subtractions non- 
perturbatively [13]. We propose to do this by 
imposing a suitable renormalization condition on 
the matrix element taken between heavy quark 
states, and to compute the corresponding sub- 
traction constant using lattice simulations. For 



example we can impose the natural condition that 

{h{p = 0)\hDlh\h{p = 0)) = (34) 

in some gauge, e.g. the Landau gauge, where the 
subtracted operator TiDgh is given by a linear 
combination of JiD^h and all lower dimensional 
ones with which it can mix 

hDlh = hD'^h-cJiv-Dh-C2hh . (35) 

For the matrix elements of the subtracted opera- 
tor between hadronic states it is the constant C2 
which is required (the term multiplying Ci can be 
eliminated by the equations of motion) , and this 
can be determined with very good precision [19] 

a^cs = -0.73 ±0.01 ±0.02 . (36) 

The matrix element of the subtracted operator in 
lattice units is then given by the difference of the 
results in equations (32) and (36). The cancella- 
tion is overwhelming and leads to very large rela- 
tive errors. Clearly a very large statistical sample 
will be required to obtain a result for the matrix 
element of the subtracted kinetic energy operator. 

9. A 

Much of the phenomenology using the HQET 

is presented in terms of the parameter A = niH — 
niQ, where nin is the mass of the heavy hght 
hadron H. In view of the discussion above the 
question arises as to what wc should take for the 
heavy quark mass itiq. The pole mass is am- 
biguous and cannot be used [2,3]. Short-distance 
definitions of the heavy quark mass can be used 
safely, but then it is not clear whether the intro- 
duction of A is useful (e.g. in that case A does 
not go to a constant as mg — > oc). 

In section 5 we have seen that even if we 
start with the bare HQET with no mass term, 
C = hD^h, the presence of power divergences and 
renormalons generates a mass term. It is there- 
fore possible to add a mass term {Smhh) to the 
action and to fix the counterterm by imposing 
a renormalization condition on the heavy quark 
propagator. Lattice simulations provide the op- 
portunity of determining the counterterm non- 
perturbatively. The procedure would be as fol- 
lows: 



Consider the theory with no mass term, C = 
hDih. Compute the heavy quark propaga- 
tor in a fixed gauge (the Landau gauge for 
example), and determine tlic cfFcctive mass 
meff(i); arriesit) = \n{S eff (t) / S eff {t + a)) . 
Results for the effective mass from simula- 
tions at 3 = 6.0,6.2 and 6.4 are presented 
in Fig. 5. 

Add a mass counterterm 6mhh to the ac- 
tion. The propagator now changes by a fac- 
tor cxp(— /.if), where /.i = l/aln(l + 6m a), 
and hence the effective mass at any time 
changes by /i. 

Impose some renormalization condition on 
the effective mass, such as 7nea{t*) = for 
some fixed t* . At small times we also have 
a perturbative contibution to the effective 
mass 

asCp Ih , 

me«{t) = T ^ ^ 

where is the one-loop contribution to the 
anomalous dimension of the heavy quark 
field (in the Landau gauge ^/j = —6). Such 
a behaviour is clearly seen in Fig. 5. 

For large values of t we are in the non- 
perturbative regime, and do not know the 
behaviour of the propagator. We know 
that there is an exponential factor from 
the renormalon singularity. If other non- 
pcrturbativc effects do not change the ef- 
fective mass at large values of f, then 
win tend to a constant as f ^ oo. Lattice 
simulations provide the opportunity to in- 
vestigate whether this is true. The results 
presented in Fig. 5 are consistent with the 
assumption that the effective mass tends to 
a constant at large times. If this is the case 
then a natural definition of the mass coun- 
terterm is to impose that the effective mass 
is zero at large times, and define A by 

K = e + n (38) 

and a subtracted pole mass rriQ by 

tuq = iriH - A , (39) 



where H is the heavy hadron [18]. A and 
rriQ defined in this way are free of power di- 
vergences and renormalon ambiguities, and 
arc independent of the mctliod used to reg- 
ulate the ultra-violet divergences. It must 
be stressed however, that this definition of 
A assumes properties of the long-distance 
behaviour of the quark propagator. 

A further check on this procedure from lat- 
tice simulations is that the result for A 
should be independent of the lattice spac- 
ing. Within our errors this appears to be 
the case. We find [19] 

A = 180 ± 35 MeV at = 6.0 (40) 
A = 220 ±55 MeV at = 6.2 (41) 

where the errors combine statistical uncer- 
tainties with those due to the calibration of 
the lattice spacing. Combining the results 
we quote 

A = 190 1 30 MeV (42) 

as our best estimate of A. 

I end by emphasising the obvious but impor- 
tant point, that when presenting a value for A, 
one must be careful in defining precisely what is 
meant by this parameter. In particular it is not 
sufficient to say that the quark mass is the pole 
mass. 

10. Conclusions 

In this talk I have attempted to summarize the 
considerable progress which has been made re- 
cently towards the theoretical understanding of 
higher twist (or dimension) corrections to phys- 
ical processes. It can now be expected that a 
sound phenomenology of these corrections can be 
developed. A central step in this is the com- 
plete definition of liiglier dimensional operators 
(we have seen that it is not sufficient to say that 
they are defined in the MS scheme). The most 
important points have been summarized in the 
overview presented in section 2, which should be 
considered as the main part of the conclusions. 
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Figure 5. The Heavy Quark Effective Masses from simulations by the APE collaboration at 3 different 
values of (3. In each case the size of the lattice is 24^ x 40. 



I conclude with a brief mention of some re- 
lated questions. In rcfs.[7-9] the question of sub- 
asymptotic corrections to the Drell-Yan process 
is studied. In refs.[7,8] it is argued that the pres- 
ence of the leading i.r.renormalon implies that the 
leading power corrections are of 0(l/(5), where Q 
is the invariant mass of the lepton pair. In ref.[9] 
it is shown that the coefRcient of this renormalon 
is zero in the large Nf limit. It is a very important 
question whether this is true beyond this approx- 
imation. In any case, from studies of the rcnor- 
malons, it is agreed that there are 1/Q corrections 
to event-shape variable in e+e" collisions. 

It is an interesting question whether there are 
renormalon ambiguities in lattice perturbation 
theory for quantities which diverge like inverse 
powers of the lattice spacing. A simple example 
is the critical mass, merit (for which the mass of 
the pion is zero) in simulations of lattice QCD 
with Wilson fermions. In the large Nf limit the 
i.r.renormalon at it = 1/2 is absent, implying that 
there is no ambiguity in merit of O(Aqcd) [21]. It 
will be interesting to establish whether this result 
survives in higher orders in the Nf expansion. 

In the large Nf limit the coefRcient of the i.r. 
renormalon at u = 1 in the pole mass vanishes 
[3], and equivalently the u. v. renormalon at u = 1 
in the kinetic energy operator also has zero co- 
efficient. In rcf.[22] this was called the invisible 
renormalon. So far it is not understood whether 
there is some symmetry which would allow one 
to conclude that this is a general result rather 



than an accident of the large Nf limit. In the 
lattice formulation of the HQET the correspond- 
ing quadratic divergence is present [13], whereas 
using Pauli-Villars or momentum-flow regulariza- 
tion the coefficient of the quadratic divergence at 
one-loop order is zero [22]. 

We have seen that the cancellation of renor- 
malons in Wilson coefRcient functions occurs be- 
tween terms which in Rnite orders of perturbation 
theory are of different orders in the expansion in 
rnq or other asymptotic variable. The higher di- 
mension operators must be defined completely. Ji 
has revisited the extraction of the gluon conden- 
sate from the lattice measurements of the plaque- 
tte V [23] 
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(43) 



Clearly before we can determine the gluon con- 
densate from the lattice measurement we have 
to subtract the perturbative series. This has 
to be done carefully as the series has i.r. renor- 
malon singularities. It is remarkable that, using 
the Langevin stochastic formulation of the lattice 
Yang-Mills theory, the Parma group were able 
to determine the perturbative coefficients up to 
eight-loop order [24]! 

Ec„ 2 1.218 2.960 
— = - -I 1 1- 
/3" p p'^ 13^ 
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9.28 34 135 563 2488 

It seems likely that these techniques will be very 
important in other applications where high or- 
ders of perturbation theory are needed to control 
the cancellation of renormalon ambiguities. Ji's 
attempts to define a gluon condensate by sub- 
tracting various resummed versions of the series 
in eq.(44) led to results which were a factor of 
more than 5 larger than values used in QCD sum 
rule applications [23]. It remains an interesting 
question, from both the theoretical and numeri- 
cal point of view, how to connect lattice deter- 
minations of the condensates with those used in 
sum rules? Ref.[25] contains some discussion con- 
cerning the definition of condensates in sum-rule 
applications. 
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